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Kennedy Space Center is home to two liquid hydrogen storage tanks, one at 
each launch pad of Launch Complex 39. The liquid hydrogen storage tank at 
Launch Pad B has a significantly higher boil off rate that the liquid hydrogen 
storage tank at Launch Pad A. This research looks at various calculations 
concerning the at Launch Pad B in an attempt to develop a solution to the 
excess boil off rate. We will look at Perlite levels inside the tank, Boil off rates, 
conductive heat transfer, and radiant heat transfer through the tank. As a 
conclusion to the research, we will model the effects of placing an external 
insulation to the tank in order to reduce the boil off rate and increase the 
economic efficiency of the liquid hydrogen storage tanks. 
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The Liquid Hydrogen storage tank at Launch Pad B has a significantly 
higher boil off rate than the Liquid Hydrogen storage tank at Launch Pad A. 
Both tanks were built in 1965 by the Chicago Bridge and Iron Company of 
Salt Lake City, Utah 1 ; however, the tank at Pad B has consistently had a boil 
off rate more than double that of the tank at Pad A. The tanks are constructed 
to keep the liquid hydrogen extremely cold; consequently, there is insulation 
between the inner tank and outer tank. The insulation that was used for 
the construction of the tanks is Perlite, which has a powdered texture. The 
outside of the tank at Pad B has also developed a cold spot where we believe 
there may be a Perlite void. We will look at several calculations in order to 
better understand the make up of the tank, its dimensions, and to explain 
why there is such a significant difference in the boil off rates between the two 
tanks. We believe that refurbishing the tank is not the only solution to the 
existing boil off problem; consequently, this research will attempt determine if 
coating the cold spot on the outside of Tank B with Spray on Foam Insulation 
will reduce the boil off rate and increase the economic efficiency of the tank. 

The tanks are constructed so that there is an inner tank and an outer 
tank. The inner tank, made of stainless steel, sits suspended on supports 
inside of the outer tank, made of carbon steel. Between these two tanks, the 
powdered insulation, Perlite, is pumped in and then a vacuum is created. On 
the outside of Tank B, there is a cold spot which tells us that there is heat 
escaping the tank from that area at a higher rate than at any other area 
on the tank. It is our belief that there exists a Perlite void at that location 
between the inner tank and the outer tank and that a possible solution is to 
use the Spray on Foam Insulation to lower the temperature on the outside 
of the tank thus reducing the temperature flow inside the tank leading to a 
lowered boil off rate. 

As a beginning to this research, we will develop a formula that will 
allow us to convert the liquid volume readings to liquid level inside the tank 
in meters. This will allow us to better understand how the boil off rate relates 
to the liquid levels in the tank. This calculation applies to both tanks and 
therefore can be used for further analysis of related topics. The following 
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equation allows us to go from volume to the height in the tank: 

irh 3 {r ~^) = V 

Where: 

h = The height from the top of the liquid to the top of the tank in meters. 

r = The Radius of the tank. 

See Appendix A for the derivation of this formula. We can now fill in the value 
of the radius 2 and use the resulting formula to solve for various heights given 
specific volumes inside the tank with a mathematical computer program: 

7t/i 3 (9.35 - ^) = V 

After the tanks were constructed, the inner tank cooled as the liquid 
hydrogen was added to the tank. During this cooling process, the inner tank 
contracted which caused a decrease in the volume of the inner tank. The 
Perlite between the tanks dropped during this process relative to the change 
in volume. In order to determine the amount the Perlite may drop due to 
the change in volume of the inner tank during cooling, it is necessary to 
calculate the change in volume of the inner tank due to thermal expansion. 
One assumption that we make with this calculation is that the Perlite acts as 
a liquid between the tanks. Given this assumption, the value that we get will 
be a slight overestimate. First, we start with the volume change of the inner 
tank. The following equation is the thermal expansion equation for stainless 
steel 3 : 

— ~ L293 = (a + bT + cT 2 + dT 3 + eT 4 ) • 10“ 5 
T293 

Where: 

a = -2.9546 * 10 2 
b = -4.0518* KT 1 

2 Berg, Mark. Heat Transfer Analysis of the LC-39 Liquid Hydrogen Storage Tanks. Powerpoint Presen- 
tation, pg. 21 

3 Cryogenic Material Properties Database ; E.D. Marquardt, J.P.Le, and Ray Radebaugh; National Insti- 
tute of Standards and Technology; June 2000 
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c = 9.4010 * 1(T 3 
d = -2.10989 * 10- 5 
e= 1.8780 * 1CT 8 
T = Temperature, (K) 

L = Length at a given temperature in meters. 

The tank cools down from room temperature (293K) to 20K as it is filled 
with the liquid Hydrogen. Therefore, the thermal expansion formula yields: 

Lt ~ -£>293 0 nr . , n _ 3 

— = —3.00 * 10 m. 

L 293 

The radius of the inner tank at room temperature 4 is given as: L 293 = 9.373m. 
Using the previous result to solve for Z/ 2 o yields a new radius of 9.345m. 

The volume of the inner tank at room temperature ( Vo i ) is: 

4 Q 

V = -7r r 3 
3 

V 0i = ^7r(9.373m) 3 
^ = 3449.25m 3 

After the tank is cooled to 20K, the new volume of the tank ( V\ t )is: 

Vii = ^7r(9.345m) 3 
V u = 3418.43m 3 

Consequently, the change in volume due to cooling of the tank ( AV ): 

AU = U 0i - V u 

4 Berg, Mark. Heat Transfer Analysis of the LC-39 Liquid Hydrogen Storage Tanks. Powerpoint Presen- 
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Ay = 30.82m 3 

AF can also be calculated using the volume between the tanks: 

Vo. = ^7r(10.65m) 3 
y 0o = 5059.85m 3 

Where Vq o is equal to the volume of the outer tank. 5 Then the volume between 
the tanks at 293K is equal to: 


v 0b = v 0 . - y 0i 

y 0i = 1610.6m 3 

and 

Vi b = Vo. - V u 
V h = 1641.42 m 3 

Where Vo b is equal to the initial volume between the tanks before cooling and 
V u is equal to the volume between the tanks after cooling. 

So the change in volume between the tanks ( Ay ) is: 

Ay = v u - v 0b 

AV = 30.82m 3 

Which is the same Ay as the previous calculations. 

Now that we know how much the volume changed due to the cooling 
of the inner tank, we can use the same formula that we used for the volume 
to liquid level conversion to determine if the perlite dropped below the top 
of the inner tank. It should be noted that this formula only works for this 
calculation if the Perlite drops an amount less than the distance between the 

5 Berg, Mark. Heat Transfer Analysis of the LC-39 Liquid Hydrogen Storage Tanks . Powerpoint Presen- 
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tanks. This calculation also applies to the tanks at both Pad A and Pad B. 
Beginning with the conversion equation: 

*h Z (r ~^) = V 

we can use the dimensions for the outer tank and volume change of the inner 
tank to determine if the Perlite drops below the top of the inner tank. We let 
r =10. 64m, the radius 6 of the outer tank, and we know the change in Volume 
is 30.82m 3 . This yields the following equation: 

7r/i 3 (10.64m — — ) = 30.82m 3 

Now, using mathematical computer software such as Mathematica, we solve 
for h and find that if the Perlite acts as a liquid, it would drop an amount 
equal to 1.0m. The distance between the two tanks is equal to 1.3m; therefore, 
the top of the inner tank was not exposed due to the contraction of the inner 
tank. Since we found that the Perlite dropped less than the distance between 
the tanks, using this method works for this case. 

An essential assumption made for this research is that there is a Perlite 
void between the tanks at the area of the cold spot. In order to demonstrate 
the validity of this assumption, we looked at the excess boil off rates of Pad B 
when compared to Pad A. Building on the research of Mark Berg, we take the 
excess boil off rate and compare it to the formulas for conductive heat transfer 
and radiant heat transfer. With this data, we can produce valid evidence that 
supports the claim that the cause of the higher boil off rated is due to radiant 
heat transfer. First, we tried to produce a constant, k, that would explain the 
heat transfer for both cases. However, we found that the constant does not fit 
the model if there is solely conductive heat transfer through the tank. This 
comparison explains why the darker grey tank had a much higher boil off rate 
than the white tank. The following chart was produced that demonstrates 
the difference in the conductive heat transfer and the radiant heat transfer: 


6 Berg, Mark. Heat Transfer Analysis of the LC-39 Liquid Hydrogen Storage Tanks. Powerpoint Presen- 
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Given these temperatures, it is not possible to produce a constant 
that explains the higher boil off rate for the tank at Pad B using only the 
conductive heat transfer. 

Further evidence that there is a Perlite void at the Pad B tank can be 
seen in inferred images. Appendix B is a thermal image of the LH2 tank at 
Pad B. This image can be compared to Appendix C which is an inferred image 
taken of scale models of the liquid hydrogen tanks where Perlite voids have 
been created. They appear to have similar qualities the lead us to believe 
there is a Perlite void at Pad B. 

In order to demonstrate the differences between the insulation provided 
by Perlite and that of a Blackbody, we will first exam the effects of placing 
SOFI on the outside of the tank if there exists Perlite between the outer and 
inner tanks. We will later compare this result to the effects of placing SOFI on 
the tank where there exists a void. By calculating the amount of heat transfer 
through the Perlite, we can demonstrate how adding varying thicknesses of 
SOFI to the tank will affect the overall heat transfer through the tank. We 
will use the following simplified model of the tank with Perlite: 
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In this model, the outer insulation is Spray On Foam Insulation (SOFI). 
The temperature of the middle of the two insulations varies with the thickness 
of the SOFI. Therefore, the thickness of the SOFI is used as the variable L. 
Where: 

T 0 = temperature on the outside of the SOFI (300K) 

T m = temperature(K) between the two insulations, variable temperature 
Ti = The temperature on the inside of the inner tank(20.28K) 

X 2 — xi = The thickness(m) of the Perlite between the tanks (1.308m) 

X 3 — X 2 = L = The thickness(m) of the SOFI, variable distance 

The equation 7 for the energy transfer for a homogeneous substance is given 
as: 

q — —kAVT ■ n 

Where: 

q = energy transfer rate 

VT-n = temperature gradient in the direction normal to the area, A 
k = thermal conductivity of insulation 

7 Pitts, Donald. Heat Transfer . McGraw-Hill Book Company. 1977. pg. 1 
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Since this is a one dimensional case, VT-n becomes Therefore, assuming 
a linear model for the heat transfer allows the following substitution: 

AT _ T 2 - Ti 
Ax X2 — X\ 


This produces the new equation: 


q = 


—kA( 


T2-T1 

X2 ~ Xi 


) 


With the simplified model, the energy transfer rate from the outside of 
the SOFI to the middle and the energy transfer rate from the middle to the 
inner side of the Perlite, must be equal. Consequently, in order to solve for 
the varying temperature and thickness of SOFI, it is permissable to set the 
equations for the two types of insulation equal to each other. This yields: 

- k sofi A ( T ° ~ Tm ) = -V^(— ~ -) 

£3 - X2 X 2 ~ Xi 

Where: 

ksofi = 2.4 * 10 
kperl = 1*5 * 10 ^ — 


Filling in all the values: 


-2.4 * 10-2( 300 = _1 5 * 10 - ^ 20 28 ) 


4.29 


T = 

± vn — 


7.09 * 10~ 3 + y 


3.50 * 10- 4 + 


2 . 4 * 10- 2 


The following graph shows the thickness of SOFI ranging from 0-10m. and 
temperature ranging from 0-300K. 
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It is evident that the temperature changes slowly with the thickness of 
the SOFI when there is Perlite between the tanks. 

Now we will develop a new model, removing the Perlite from the model 
and replacing it with a blackbody state. Since our previous calculations sug- 
gest that this is a more probable explanation of the heat transfer through the 
cold spot, we can examine the effects of adding SOFI to a tank with these 
conditions. 


SOFI 



Blackbody 


For these calculations, we will use the formula for heat transfer in a 
blackbody. The derivation of this formula can be found in Appendix D. In 
order to examine the heat transfer through the tank caused by radiant heat, 
we will assume that the outer tank has an emissivity of 0.2 and the inner tank 
has an emissivity of 0.8. Using the heat transfer formula in Appendix D, the 
net emissivity of the tank is 0.19. Now, as with the previous calculations, we 
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set the two heat transfer formulas equal to each other to determine the effect 
of the SOFI on the heat transfer for the tank: 


T —T 

k S ofiA (— — — ) = 
X3-X2 


~ aTf) 

62 + 6i(l — 62) 


By filling in the values that we know, we get the following equations: 


300 — T 

2.4 * 10" 2 ( —^) = .19 * (5.67 * 10“ 8 )(T^ - 20 4 ), 

L 

Now we can solve this equation for L and produce plots that will relate the 
thickness of the SOFI on the outside of the tank to the temperature of the 
outer tank (T m ). 

2.4 * 10 12 - 8 * 10 9 T m 
“ -5.7456 * 10 8 + 3591T4 

In order to compare the effect of adding SOFI to the outside of a tank that 
has Perlite to the effect that it has on a tank with a Perlite void, we plot the 
two graphs together: 



Temperature (K) 

Perlite Void Perlite Present 

It appears evident from this plot that adding SOFI to the outside of 
the tank would reduce the temperature of the outer tank by a greater amount 
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if there exists a Perlite void. This would also lower the amount of heat that 
is transferred through the outer tank to the inner parts of the tank. The plot 
also demonstrates that SOFI would not have as much of an effect on areas of 
the tank where the perlite is in place. 

There is much more research to be conducted in this area. Using scale 
models of the LH2 tanks to test this hypothesis would be a good direction for 
future research. Another area that needs examining is a more complex anal- 
ysis of the heat transfer taking into consideration the heat transfer through 
the steel and into the outer edges of the SOFI. 
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Appendix A: 

To find an equation to go from the volume of a sphere to the liquid 
level inside the sphere, we begin with the general general equation of a circle: 

x 2 + y 2 = r 2 


From this, we solve for x 2 : 


222 
x —r—y 


Now, we can use an integral to represent the volume: 

7 r / r 2 — y 2 dy — V 

Jr—h 

This becomes: 

n(r 2 y - ^)\ r r _ h = V 

Now, filling in the limits of integration, yields: 


Simplifying: 


7r[r 3 — - — [r 2 — (r — h) — 


- h) [3r 3 - (r - h) 2 ]} = V 


(r - hf 


]] = V 


Now, we can factor out a 4 : 


7T 


-^■[2r 3 — (r — h)[3r 3 — (r 2 — 2 rh + h 2 )]] — V 
o 


Which is equal to: 


7T 


'-{2 r 6 - 2 r 6 - 2 r 2 h + rh 2 + 2r 2 h + 2rh 2 -h 3 ) = V 


Combining like terms yields: 


7T 


^(3 rh 2 -h 3 ) = V 


Factoring out the h 2 : 


irh 2 


(3 r -h) = V 


And this simplifies to our equation for converting volume to height: 

nh 2 h 
— (r -3 ) = K 


12 



Appendix B 



17.2 D C 


Appendix C: IR Photographs 
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Appendix D: 

The following is the derivation of the heat transfer equation for a black- 

body: 


Q — > 


Surface 1 



Surface 2 


El 



E2 



Blackbody 


The net energy lost from surface 1 (E\ net ) is given by the directly ra- 
diated energy minus the reabsorption of the reflected energy and the corre- 
sponding radiation from surface 2. E is the radiant energy flex per unit area 
and is given by the Stephen- Boltzmann Law E = eaT 4 . Lower case epsilon, e 
is the surface emissivity and sigma, a is the Stephen-Boltzmann constant of 
5.67 * 10 _8 ^4- Rho, p is the surface reflectivity. The relationship between p 
and e for an opaque surface is p + e = 1. 

El net = Ei—Ei€iP2 — Ei€ip2plP2 — E\€\p2P\P2P\ • • ■ — {E2£\+E2€\P2P\+E2€\P2PlP2P\ • • •) 


Rewriting the formula in summation notation yields: 

00 OO 

= [£i(i -dE PM' 1 )] - [£*1(E P 2 P 1 )} 

n= 1 77 —O 

Changing the index on the first term and simplifying the second term: 

OO OO 

= [£1(1 -€i E P? + V?)] - [JWEOwin 

n = 0 n=0 
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Simplifying the first term: 


OO OO 

= [£i(i - ei £ wM)] - [B2£i(E M n \ 

n— 0 n— 0 

Pulling out the constant p2‘. 

OO OO 

= [Ei(l - eip 2 EM n )} ~ (p 2 pi) n ] 

n= 0 n= 0 

Using the identity l+a:+a; 2 +a; 3 . . . x n = for |rr| < 1, the formula becomes: 

-)] 


= [£,(1 - £i( , 2 (-J—))] - [E 2 ei(- 1 

1 — jOojOi 1 


P2Pl 


Distributing: 


= i^i(i 


€\P2 ^ -^ 2^1 j 

1 — P2Pl 1 — P2Pl 


Obtaining a common denominator for the first term: 


[El( l f 2 *’ 1 6lP2 )] - [: 

1 — P2Pl 


E 261 
~ P2Pl 


Distributing: 


_ Ei — E 1 P 2 P 1 ~ Eie\p2 ~ E^e 1 
1 - P 2 P 1 

Substituting p\ = 1 — £1 and P 2 = 1 — € 2 : 

_ Ei — Ei( 1 — £ 2 ) (1 — Ci) — £i€i(1 — € 2 ) — E 2 G 1 
1-[(1-£ 2 )(1-£ 1 )] 


Distributing: 

Ei — Ei + Ei 6\ + £’162 — £iCi£ 2 — £161 + £ieie 2 — £ 2 e i 
1 - (1 - £1 - e 2 + €261) 
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Canceling terms yields: 


__ E\€2 — i?2 e 1 

Cl + €2 — 6162 

Simplifying the denominator: 

E\€2 — E 2 C 1 

e 2 + ei(l-e 2 ) 

Substituting E\ = eicrT 4 and E 2 = € 2 ^ 2 : 

ei<7T 1 4 C2 — e 2 crT 2 4 ei 
€2 + €l(l — ^ 2 ) 

Simplifying: 

6 1 € 2 (crT 1 4 — crTl) 
e 2 + ei(l - e 2 ) 

This formula can be found in cryogenic texts such as Flynn, Thomas; Cryogenic Engineers 
pg. 363. 
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